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ABSTRACT 

Several techniques for finding approximate solutions to 
certain classes of nonlinear differential equations are in- 
vestigated. The nonlinear systems evaluated are second order 
with quadratic and cubic terms, having driving functions and 
initial conditions. Pipes' technique is used to reduce the 
nonlinear differential equation to a system of linear dif- 
ferential equations. The Brady-Baycura technique makes use 
of nonlinear Laplace transforms to obtain the solution. The 
solutions are compared to the well known Runge-Kutta numerical 
method solution using the digital computer. Greater accuracy 
was found using the Brady-Baycura method, but the simplicity 
of Pipes' method makes it more attractive to the engineer. 
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I. 



INTRODUCTION 



Although numerical methods such as Runge-Kutta have been 
widely used to solve nonlinear differential equations by means 
of digital computers, the engineer does not always have access 
to a computer. Practical techniques to assist the engineer in 
obtaining approximate solutions are needed. 

Pipes (Ref. 1) developed a method of solving nonlinear 
differential equations by a reversion technique quite akin 
to algebraic series reversion. The technique uses Laplace 
transforms extensively in obtaining the approximate solution. 
Baycura and Brady (Refs. 2 and 3) have proposed and examined 
nonlinear Laplace transforms to obtain solutions to nonlinear 
differential equations. Tou , Doetsch , Pipes and Nowacki 
(Refs. 4, 5, 6 and 7) have proposed an iterative method using 
Laplace transforms in each iteration. Flake (Ref. 8) has pro- 
posed a cumbersome method using Volterra series and special 
transforms . 

The techniques were investigated and applied to driven 
second order systems with initial conditions. The systems 
were restricted to known stable systems, and with driving 
functions far from resonance. The method of Flake was found 
to be so impractical that it is only mentioned here for 
reference . 
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II . METHODS 



The various methods of obtaining approximate solutions to 
nonlinear differential equations were investigated. 



A. PIPES' METHOD 

Pipes has postulated a general nonlinear differential 
equation of the form 



a* v- a A v<S a ( x l + . . . a n x h = k <f>w 



(2.1.1) 



where t is the independent variable, x is the dependent vari- 
able, k is a constant, 0 (t) is a driving function, and the a^ 
are functions of the differential operator D, where 

D = d/dt (2.1.2) 



Consider the special case of equation (2.1.1) where 




d=o 




and the a^ are constant coefficients for i - 2, 3, 
Equation (2.1.1) becomes 



( + L>, 1> <ii_ + 

at ^ 




X + Q^X 2, + . , . 

+ Qn X r '= b fiO?) 



(2.1.3) 



n. 



(2.1.4) 



or using the operator notation 

a„ + W,D +k A D+ . . . + aX+ • - • 

+ Q„x" = b 0Cf) 



Following Pipes, assume a series solution of the form 

xct) = + Acok z +Amk 3 + . . . ( 2 .i. 6 ) 



6 



By substituting equation (2.1.6) into (2.1.1) and equating 
coefficients of equal powers of k, the various A^t) terms 



may be found. They are 

(2.1.7) 

/\ 2 W- - a fc A*U:)/ a . (2.1.8) 

A 5 (0 - ‘ tMO ± (2.1.9) 

A,W - - OA)[ol.(A> ' 2 ^ A 1 , to A,« + a., A* w 1 (2.1.10) 



f\ 5 ^>- ~ Aa* z>c\ 2 (j\(X) Aj_oo 4 

+ 4a 4 ^,wA z W 

Pipes (Ref. 1) lists additional terms. However, since a 
second or third approximation to the solution will generally 
be all that is desired, only the first five terms are given 
here. 

A first order approximation can be found by solving the 
linear differential equation (2.1.7). Since (2.1.7) is a 
linear differential equation, it can be solved by the usual 
techniques of integration or Laplace transforms, applying 
initial conditions. The remaining A^(t) are solved in like 
manner, except the initial conditions are zero. 

B . BRADY- BAYCURA METHOD 

Brady and Baycura have shown that for exponential 
functions and for small values of time the following approxi- 
mate expression holds 

£[xVi]~ (2.2.1) 
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Since the Laplace transform of 0(t) is known, all trans- 
forms form the terms of equation (2.1.3) are known a..d thus 
(2.1.3) can be solved for X(s). 

C. TOU-DOETSCH-PIPES-NOWACKI METHOD 

This method probably originated with Pipes, but was 
arrived at independently by the others. 

Assume the nonlinear system differential equation of the 

form 




Taking the Laplace transforms of both sides of equation 
(2.3.1) yields 

r> 

X /^> Qo} ^ 0: ^ y . .'jj ~ GyCv) (2.3.2) 

where Q(s) is the set of initial conditions. 

In the first approximation assume the nonlinear term 
then 

y W . 4 ^ ‘ (2.3.3) 

The inverse Laplace transform of Y^(s) yields the time 

domain solution of the first approximation. 

The nonlinear term is included in the second approximation. 

y 4 ‘ -0 

• — (2.3.4) 

k - O 
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. . . ) is a function of the first approxi- 



where f(t, y^, y^, 
mation of y^ (t) . 

Successive approximations are found by continuing the 
same procedure. 
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III. APPLICATIONS 



Nonlinear systems investigated were limited to second 
order with initial conditions and extended to systems with 
forcing functions. 

A. NONLINEAR SPRING 

Consider a simple spring-mass system described by 

m -v Fw = o (3.1.1) 

The system is undamped, with an initial displacement of 
x(0) = x q/ and initial velocity set to zero. F(x) is the 
spring force, and for this example, let its characteristic 
be 

= kx + (3.1.2) 

Letting all constants be unity, the system equation is 
nonlinear 

di. + X -v- X z = O (3.1.3) 



» o 



(3.1.4) 



1. Solution by Pipes Reversion Method 

Using the method outlined in Chapter II, Section A, 
equation (3.1.4) can be rewritten as 
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(D 2 + 0 x 4- X z = o 



(3.1.5) 



The a. coefficients are seen to be 
i 

Q,. - i l 

k - 1 

<2(0 ~ O 



(3.1.6) 



The individual coefficients are then computed using the 
formulas (2.1.8), (2.1.9) and (2.1.10): 

-- sAo/g, . 



Cd 7 vi)a,= ° 



(3.1.7) 



Taking the Laplace transform of equation (3.1.3) and 
impressing initial conditions yields: 




sy. 

.s 1 + l 



(3.1.8) 



This well known transform gives 

A, CO - 



(3.1.9) 



As expected, the system is oscillatory. 

The next term A 2 (t) i s f° und from equation (2.1.8) 




-_xl_Co£jt 

V? 4- | 



(3.1.10) 



A 2 (t) and all subsequent A ± (t) are determined by the con- 
dition that they have zero initial conditions. 
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Using the identity 



Cos 7 t - -L(j +OsZ.-t) (3 

and taking Laplace transforms of both sides of equation 
(3.1.10) 



Both terms can be expanded into partial fractions: 



A,0> -T \i - s vq 



'h-> 

^ + I 



+ 4 



1 



(3. 



Taking inverse Laplace transforms: 

\(& = ~ -i-coszt^ (3. 

A^(t) is found in a similar manner, and the solution 
given by 

4 \ A 3 00k\ . . , 

x(fc)«X£.oU " V l - 1 c*s\ - ^ Cj^ X’t*! 

2_ L 3 J 

^ X » W ^.S\rA- ^ 4 i!L Co^4-)-Ct»S2rt+-i- tos'S'll 4 . , v ( 3 . 
3 I 1 a-s 3 it. -J 



In equation (3.1.15) it is observed that the first 
secular 1 term occurs in the third approximation. 



1 Secular terms are those involving t sin t, t cos t, t 
. . . , etc . 



1.11) 



1 . 12 ) 



1.13) 



1.14) 

is 



1 . 6 ) 



1.15) 



sin t , 
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2. Solution by Brady-Baycura Transform Method 



Taking the Laplace transform of equation (3.1.3) and 
applying initial conditions gives 

-SOC* + Xts5 ^ - O (3.1.16) 

In the usual manner, equation (3.1.16) is solved as a 
quadratic. After rearranging terms 



Yes} - 



sH \ 




(3.1.17) 



Neglecting the negative radical as physically impossible, the 
positive radical can be expanded as a binomial 



Xca - 









(3.1.18) 



The first term is just x Q cos t; however^ Lcs^O^J 

is not one of the common transforms. Reference 10 lists the 
following transform: 



z 












(3.1.19) 



where Jn-3/2(w Q t) is the half order Bessel function. 

Further it is known from the theory of Laplace trans- 
forms that 



13 



\?=c> 



v\~\~ Vi 

<j* -b 



n= . 



(3.1.20) 



Since f (0) = 0 for k = 0, 1, 2, . 
can be found : 

letting £>„= \ 

J y « - (vt)H ^ -“**) 



, n-l, the second term 



(3.1.21) 



and after carrying out the differentiations of equation 
(3.1.20), the second term of equation (3.1.18) becomes 






*»■ -v 

s L 






(3.1.22) 



Repeated application of the method will yield the 
third term. 

The approximate solution to three terms is the truncated 
series 

* D Oovt 

+ (3.1.23) 
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3 . Solution by Tou-Doetsch-Pipes-Nowacki Method 

Putting equation (3.1.3) in the form of (2.3.1) 



0, y * (XX \ Y? = O 

«- o 



(3. 



where 



C ' t> - ' 
a, - o 



3 (A)=o 



J 



(3.: 



Taking Laplace transforms of both sides 






(3.: 



Q- * -V Z 1 “ 



o 



(3 . J 



Solving for X 1 (s) as in equation (2.3.2) 



X.<3 - 



Qd^ 

t ^ 



f> X 4 \ 



(3.] 



Then taking inverse Laplace transforms yields the first 
approximate solution: 



X b CLD^ 



(3.] 



1.24) 



. 25) 



. 26) 



. 27) 



.28) 



.29) 
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The second approximation is found by applying equation 



(2.3.3) 






V-*> 



(3.1. 30) 




•S^V\ 






\x~ 



(3.1.31) 



Using the identity of equation (3.1.11), equation 
(3.1.31) becomes 












^ V4 - J 



(3.1.32) 



Taking the inverse Laplace transform gives the second approxi- 
mation as 

X, ©* - £ L _ Z. . i 2-t] 

~i__ 2> I 



(3.1.33) 



Repeating the method gives x^(t) and the result is the same 
as equation (3.1.15). 

While this method appears to be different than Pipes 
method of Chapter II, Section A, consideration of the manipu- 
lation of the solution shows the two methods to be equivalent. 
Since the method of Chapter II, Section A is less cumbersom, 
it will be used hereafter. 
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4. Comparing the Methods 



Figures 1 and 2 are phase portraits of the solution 
as given by the IBM 360 computer. Runge-Kutta integration 
technique was used to obtain the solution. From these plots, 
it is clear that the system is oscillatory for x £ 0.5 units 

o — 

and unstable otherwise. 

Figures 3, 4 and 5 compare the solutions by Pipes 
method and the Brady-Baycura technique for various initial 
conditions . 

The graphs demonstrate the relative accuracy of the 
approximations. The Brady-Baycura technic gives a better 
approximation in each instance; however. Pipes method is as 
accurate for small initial conditions. The difference is 
caused by the location of the secular terms in the solution. 

For both methods, however, the approximations are of little 
value for times greater than one unit. 

No attempt has been made to suppress the secular terms 
which cause the approximate solution to diverge. It is gen- 
erally known that subharmonic resonance terms (submultiples 
of the natural frequency) exist in these nonlinear systems; 
they do not occur in this solution. The higher frequency jump 
resonance terms do occur. 

B. NONLINEAR SPRING WITH DAMPENING 

The system of example 3.1 was modified to include frictional 
dampening. The resulting system equation takes on a velocity 
term. Letting all coefficients be unity, except the velocity 
term, the system equation becomes 
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FIGURE 1. Phase portrait of x + x + x - 0, x Q 1.0 



18 




FIGURE 2. 



Phase portrait of x + x + 



2 



x 



0 , 



x 



o 



0.5 
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3. Time solution of x + x + x 0 > x o 
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4. Time solution ofx+x+x “ 0, x q - 0.5 
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FIGURE 5. 



Time solution of x + x + x 



0 , x _ 0.25 
o 



22 



TABLE I 



Time 


X o = 1 

Runqe-Kutta Brady 


Pipes 


0.000 


1.00000 


1.00000 


1.00000 


0.050 


0.99750 


0.99750 


0.99854 


0.100 


0.99002 


0.99002 


0.99418 


0.150 


0.97762 


0.97762 


0.98697 


0.200 


0.96039 


0.96038 


0.96424 


0.500 


0.76476 


0.76469 


0.86483 


1.000 


0.20173 


0.22176 


0.56088 


1.500 


-0.43604 


-0.23227 


0.24276 


0.000 


0.50000 


x = 0.50 
o 

O.50OOO 


0.50000 


0.050 


0.49906 


0.49906 


0.49919 


0.100 


0.49626 


0.49626 


0.49678 


0.300 


0.46675 


0.46674 


0.47137 


0.500 


0.41002 


0.41002 


0.42252 


1.000 


0. 17967 


0.18245 


0.22407 


1.5000 


-0. 10512 


-0.075097 


-0.02327 


0.000 


0. 25000 


x^ = 0.25 
o 

0.25000 


0.25000 


0.050 


0.24961 


0. 24961 


0.24963 


0.1000 ' 


0.24844 


0.24844 


0.24850 


0.500 


0.21213 


0.21213 


0.21369 


1.000 


0.11163 


0.11214 


0.11715 
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(3.2.1) 




XCrt = 

%< 3 >) = O 



(3.2.2) 



The investigation included solving the system equation by 
Pipes and Brady-Baycura techniques. The final solution is 
found by letting the coefficient C approach zero in both 
methods . 

1. Solution by Pipes Method 



By letting C = 2 jb , the algebra is simplified and the same 
result will obtain by letting 2 go to zero. 

Proceeding with the solution = 0 / and upon 
substituting and rearranging 



Using operator notation, the aVs are 




(3.2.3) 




(3.2.4) 



Taking Laplace transforms and impressing initial conditions 



equation (3.2.4) becomes 
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Cr? - A- CS + \) A (>A 



(3.2.5) 



or 



A.w ' 



x.(mQ 



-i^7 X 5 V'VJK i ) 



(3.2.6) 



Reference 11, formula 119, gives the inverse trans- 
form of equation (3.2.6) and 

A,« ' "p M ■ 

s — , (3.2.7) 



letting 2^ go to zero 



AC-O- K. W (,-^j -Vexp^ 'j 



(3.2.8) 



A to* 



(3.2.8a) 



The solution for A^, after a modicum of algebra 



becomes 



A/o - ~ -A ^ cbss^ - -l ct>sz4.”J 



(3.2.9) 



and substituting ^(t) into equation (2.1.5), (t) becomes 



- \ Cc-sA. - ^ Cos ^ - r 






(3.2.10) 
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or 



A 



t \ ^ 

& = s 



U+Z/S’D 



4| 



U1 A— 

3. 1= S <•= 



(3.2.11) 



Taking Laplace transforms, equation (3.2.11) becomes 



\ 



Vc ^ 



^ l + 2 ( 4^' v 0 (A 1 ' '•*■ 01> Z + 2^3 w i) 

Vs ^ 



(3.2.12) 



Vt S 



(^-v-V^+Z /5s» + 0 + 2/2^ lV 



All inverse transforms of the terms of equation (3.2.12), 
except the second, can be found in Ref. 10. The second term 
is evaluated by convolution as follows: 

Let 



*,<a= -?T7 



(3.2.13) 



2^0 “ 



s ’--Vz,iS+ \ 
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(3.2.14) 



then z^(t) ~ cos t and z^(t) ~ sin t when 2 y6 goes to zero. 
Convolving z^ and 



6 

-t 



£>\o "t. 



(3. 2.15) 



The approximate solution to three terms by Pipes 
technique is 
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k (■t') x — < y ( 



L. 




* J 




-3 v + .. . (3.2.16) 



4-E. •s tfa. 



It is observed that the first secular term occurs in 
the third approximation. 

2. Solution by Brady-Baycura Method 

The solution proceeds by taking the Laplace trans- 
form of equation (3.2.1), and impressing the initial con- 
ditions. The nonlinear Laplace transform is used for the 
2 

x (t) term. 



equation (3.2.17) can be made into a series by using the 
formulas of Ref. 9. 




(3.2.17) 



By reversion of series technique, the left side of 



Let equation (3.2.17) have the form 




(3.2.18) 



where 




(3.2. 19) 
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then by reversion of the series 



xuV- 15^ 




(3. 2. 20) 



where 






(3. 2. 21) 




V 



(3.2. 22) 




s 

-+C-<b>-V l^” 2. 




2 b 7 

(s'-t -V t 



or letting C ~ 




X UV 

+ \ 
-V 



(S-\Zfti> + 0 3 
(z?-*Z/3>!> \ O s 



(3. 2. 23) 



(3.2.24) 



(3.2. 25) 



The first term of equation (3.2.25) is the same as 
equation (3.2.6), and is just x q cos t, as in the Pipes 
solution. 

Evaluation of the second term must be found by con- 
volution : 



Let 
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s> 

5?- -VZf3 + * 



(3.2.26) 



*?•* i(i + i 



(3.2.27) 



Taking inverse transforms, and letting 2^ go to zero, 

do&t (3.2.28) 



=• (3.2.29) 



Then the inverse of the second term is the con- 
volution 



Zcw = *-2b t C^ £ ^c-o 



(3.2. 30) 



Again after a modicum of calculus, algebra and trigonometric 
substitutions, equation (3.2.30) reduces to just 







-V- "t 



(3.2.31) 



and to two terms, the Brady-Baycura technique yields 

KC-0 ~ '^ & <Ld£ =^ - ^ co^\ v 4.'^ ^Tj v (3.2.32) 
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The technique has lost its desirability even at this 
point, since the inverse transforms are so difficult and time 
consuming to evaluate. 

3 . Comparing the Solutions 

The phase portraits given in Figures 6, 7a and 7b 
show the apparent warping of the phase plane due to the 
inclusion of the velocity term. The instability of the 
system is seen to depend on both the value of C and x . 

The Runge-Kutta solution is compared with the approxi- 
mate solutions of equations (3.2.16) and (3.2.32) in Figures 
8, 9 and 10. Table II lists selected values of the Runge- 
Kutta solution and the approximations. 

It is clear from the table that neither approximations 
are good even at 0.5 time units. It is further obvious that 
for values of the coefficient CtZ.1.0, the system equation 
(3.2.1) will approach equation (3.1.3). Under the condition 
of C small, the techniques might have some value in approxi- 
mating the transient state for small times (milliseconds) . 

C. NONLINEAR SPRING WITH CUBIC TERM 

Consider again equation (3.1.1), but let 

Y'LA- tx-vW 1 (3.3.1) 

Also let all coefficients be unity and the system have the 
same initial conditions as before. The nonlinear expression 
describing the system is 
£jL v X -v X* ~ ^ 

cK v (3.3.2) 
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2 

FIGURE 6. Phase portrait of x + Cx + x + x “ 0, x q ~ 1.0 
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2 — — 

FIGURE 7a. Phase portrait of x + 0.5x + x + x = 0, x Q “ 1 
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FIGURE 7b. Phase portrait of x + 0.5x + x + x^ “ 0, x q “ 1.2 
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FIGURE 8. Time solution to x + x + x 



+ x 2 “ 0, x “ 1.0 
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FIGURE 



x(0 




9. Time solution of x + x + 



x + x 2 " 0, x. 



= 0 . 
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X(t) 




FIGURE 10. 



2 - 

Time solution ofx+x + x + x 0, 



0.25 
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TABLE II 







X 

0 

II 

H * 

• 

O 




Time 


Runqe-Kutta Brady 


Pipes 


0.000 


1.00000 


1.00000 


1.00000 


0.050 


0.99754 


0.99719 


0.99667 


0.100 


0.99035 


0.98877 


0.98669 


0.150 


0.97870 


0.97478 


0.97014 


0. 200 


0.96290 


0.95530 


0.94711 






x = 0.50 
o 




0.000 


0.50000 


0.50000 


0.50000 


0.050 


0.49908 


0.49898 


0.49896 


0.100 


0.49638 


0.49594 


0.49584 


0. 150 


0.49200 


0.49089 


0.49066 


0. 200 


0.48604 


0.48384 


0.48344 






x ~ 0.25 
o 




0.000 


0.25000 


0.25000 


0.25000 


0.050 


0. 24962 


0. 24959 


0.24960 


0.100 


0.28849 


0.24836 


0.24839 


0.150 


0.24666 


0.24632 


0.24638 


0.200 


0. 24417 


0.24347 


0.24357 
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A 

XCcJi- O 



(3.3.3) 



1. Solution by Pipes Technique 



Using operator notation, the a^ coefficients are 



identified as 



a, - l^-V \ 

CX V -- D 
' 

V - ' 

QSCO- o 



(3.3.4) 



The A^(t) coefficients are found next. From equation (2.1.7) 



A,co = <3*00/ a, 



(2.1.7) 



or 

CD •+v^A >> (+') — T> (3.3.5) 

Taking Laplace transforms and impressing the initial 
conditions on A^(t) the solution is immediately seen to be 
the same as that of the quadratic case 

(3.3.6) 

Substitution of a£ into equation (2.1.4), gives coefficient 
A^ as zero. A^ is found from equation (2.1.4): 
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(2.1.8) 



A^v- 






(3.3.7) 



and after substituting for cos t, equation (3.3.7) becomes 



A CO - - ^ \ c.os,~S -C 

L TA X -V ^ J 



(3.3.8) 



Taking Laplace transforms of both sides 



A 3 (0) = 



- 



4- 1_ 



(3.3.9) 



The inverse transforms for the factors of equation (3.3.9) 
are found in reference 11, and the solution of A^(t) is 



A,« = - 0 \ £>\cA -V ^_Cx>A- 

3 -o [. A A- i 



(3.3.10) 



Since all even a^ and thus even A.^ are zero, A^ is 
zero. The coefficient A^ , after substitution, reduces to 



A s w- 



(3.3.11) 
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Substituting for the coefficients and manipulating trigono- 
metric identities, 

A-KV _5^L_ \ ^ ° "2>V -V 3-t Sx*A + -^CLc>«A 

* L 

— 1-C£>s> - Jl 



1 



(3.3.12) 



and taking Laplace transforms 






"bi. 









C=»s 



^') 3 



A -s 



- 



fr* 1 



(3.3.13) 



The last three terms have readily available inverse trans- 
forms; however, the first and second terms require special 
techniques . 

First term is evaluated by convolution: 

Let 



-2-, tO 



-A.CO 



r[ 

JCt 



S 

S'-VS 



] 



Cci-b.S-t 



\ 




3-t 



(3.3.13a) 



(3.3.13b) 



then 



(3.3.13c) 






“ -3_ totv - 2-^a CxjV 2>\ - 3 4 £>\r^ 
3^ °\t> 2 - 4 - 



(3.3.13d) 



The evaluation of the second term is found in Appendix A: 






(3. 3. 13e) 



The remaining terms of (3.3.13) are found in reference 10. 
The result for A^(t) becomes 

^ r_\ f 2.3 CCft - Z. <403334, V Cjb^S-0 

^ 32, [> V r 

-Vt(i>vr\4.- — 3\<->34^ ~ — 

^ ^ ^ _\ (3.3.14) 



The approximate solution by Pipes method is the truncated 
series 







-4 [l_ 

32, P 10 



2 3 Cr>s3: - v 







■+ - -S_ 4 . . . 



(3.3.15) 
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It is seen from equation (3.3.15) that the first 



secular term occurs in the second approximation. 

2. Solution by Brady-Baycura Method 

Upon taking Laplace transforms and using the non- 
linear expression of equation (2.2.2) and impressing the 
initial conditions, equation (3.3.2) becomes 






(3.3.16) 



Considering equation (3.3.16) to be a series, the method of 
reversion outlined in Chapter III, Section B, paragraph 2 
can be used here also to obtain all the terms: 






s'- A- \ 



(3.3.17a) 



- o 








The solution in terms of X(s) is 






a^-vi 






* «y 






(3.3.17b) 



(3.3.17c) 



(3.3.18) 



The inverse transform of the first term is just x q cos t. 

The second term must be solved using the technique outlined 
in Brady (Ref. 2) and in the example of Chapter III, Section 
A, paragraph 2. The method is straightforward, giving the 
second term as 
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(3.3.19) 



— V \£-t :s.wV. - -V^-Ccosd. i 

L J 



The second term (3.3.19) yields the higher order secular 
2 3 

terms t cos t and t sin t which will be small contributions 
for t li- 
lt is also known from equation (3.1.19) that 




Z' K PW • 1 



(3.3.20) 



By observing the result in second term, one sees that 

g 

even higher order (t ) secular terms will result from the 

third term of equation (3.3.18), and further that the 

coefficient for this third term is small (less than 
-4 

. 2 x 10 ); therefore, little is gained by solving equation 

(3.3.20) . 

The approximate solution by the Brady-Baycura 
technique is 

(3.3.21) 

3. Comparison of Solutions 

Figures 11 and 12 show the phase plane solution for 
equation (3.3.2). Here the solution is seen to be stable for 
all x q . The system always remains in a limit cycle. 
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It is also observed from Figures 13, 14 and 15 that 
for larger initial conditions, the Brady-Baycura method gives 
a better approximation to the solution for early time inter- 
vals. For the case of x q small, Pipe's method gives better 
results. This is caused by the location and number of 
secular terms appearing in the solution, which causes it 
to diverge. 

Selected values for the graphs of Figures 13, 14 and 
15 are given in Table III. 

D. NONLINEAR SPRING WITH DRIVING FUNCTION 

Consider the system of Chapter III, Section A, but now 
as a forced system. The system nonlinear differential 
equation is 

\ v\cV~ * + = \S Qjd^ uvt 

d-V- (3.4.1) 



where 

Vbd = -v (3.4.2) 

Letting all coefficients be unity, equation (3.4.1) becomes 



-i, K A 

o\A 



* Co} -- * 0 
* Cb 



(3.4.3) 



(3.4.4) 
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FIGURE 11. Phase portrait of x + x + x 



0, 



x 



o 



1.0 
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3 - 



FIGURE 12. Phase portrait of x + x + x 0, 



x =0.25 
o 
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FIGURE 13 . 



Time solution of 



x 



0, x. 



1.0 
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*(t) 




3 . 



FIGURE 14. Time solution of x + x + x 0, 



x =0.50 
o 
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XCt) 




FIGURE 15. 



Time solution of x + x 




0, x Q - 0.25 
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TABLE III 





O 






Time 


Runqe-Kutta 


Pipes 


Brady 


0.000 


1.00000 


1.00000 


1.00000 


0.050 


0.99750 


0.99768 


0.99750 


0.100 


0.99003 


0.99073 


0.99002 


0.150 


0.97767 


0.97920 


0.97761 


0 . 250 


0.93877 


0.94277 


0.93831 


0.500 


0.76880 


0.78095 


0.76278 


1.000 


0.23370 


0.25179 


0.19357 




x ~ 0.5 
o 






0.000 


0.50000 


0.50000 


0.50000 


0.050 


0.49922 


0.49922 


0.49922 


0.100 


0.49688 


0.49690 


0.49688 


0.250 


0.48064 


0.48077 


0.48063 


0.500 


0.42461 


0.42498 


0.42444 


1.000 


0.22682 


0. 22707 


0. 22681 




x = 0.25 
o 






0.000 


0.25000 


0.25000 


0.25000 


0.050 


0.24967 


0.24967 


0.24967 


0.300 


0. 23815 


0.23816 


0. 23815 


0.750 


0.17928 


0.17930 


0.17931 


1.000 


0.12946 


0.12947 


0.12966 


1.500 


0.0089567 


0.0089509 


0.011029 
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1 . Solution by Pipes Method 

The a_^ coefficients are those of equation (3.1.6), 
except here it is seen that 0 {t) = B cos^t. 

The solution proceeds as before 



Aw = <&(£/c{ x - 



1?--v \ 



(3.4.5) 



Taking Laplace transforms and impressing initial conditions. 



A ^ + 



(3.4.6) 



The terms of A^(s) have transforms in Ref. 10. 






(3.4.7) 



= W* Cosjl - 



(3.4.8) 



where 






WA-\ 



)_ - ^ 



(3.4.8a) 



The solutions to A 2 (t) and A 3 (t) are straight orward 
using the method outlined in Chapter II, Section A. Thus 

K$) - - ^ k ^Ctost- ^ 

+ 1, (3.4.9) 
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where 



V 


2 




1 


= 


^3 - ^ -V ( J_ , \ ^ 








^ tO % W-2- 












> 




6k ^'>(po~0 






■= 
















k-, ' 














J 



The solution for A^(t) is given in Appendix B. 

■ The complete solution by Pipes method is given by 

* 1 . 0 - v V v , v (3.4.10) 



2 . Solution by Brady-Baycura Technique 

Equation (3.4.3) can be transformed as in Chapter 

III, Section A, paragraph 2 using the nonlinear transform 
2 

for x (t) . After impressing the initial conditions, the 
transformed system equation becomes 



SX (.*>■) + -v ^ • ’ , (3.4.11) 

Letting equation (3.4.11) be a series, then by the method of 
algebraic series reversion 

V ' , » (3.4.12) 

s*-v^ j V J 



where the A. coefficients are 
r 
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^ + \ 



A, - 

^ 7 , = 



(s'-* 



2^ 



(3.4.13) 



Thus 

X<&> 








(3.4.14) 



The first term of the series in equation (3.4.14) 
is easily evaluated by performing the algebra indicated and 
using the transform pairs of Ref. 10. Let the first term 
be (s) , then 




-v i 






(3.4.15) 



and 



' ^6^-OA.V V — Vd£)i\ - 

\_ 6 J 



(3.4.16) 



= k.cxs^L - 



(3.4.17) 



where 
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(3.1.18) 






o 









The result is not unlike the first approximation 
obtained in the previous method by Pipes. 

Let the second term be F 2 (s) / then 



- 








(3.4.19) 



^ v 'i^V 1 



(3.4.20) 



The terms of ^ 2 ( 3 ) are all uncommon transforms. The 
technique of finding the inverse transforms is as follows: 
Since 



-C IjpWb 



(3.4.21) 



As shown in Appendix A, the derivative method of finding 
inverse transforms is extended to equation (3.4.21). 




Jl_ 

e 



SvaA: V "V" (UiVt ^ 



(3.1.21a) 



and the inverse transform of the first term of equation 
(3.4.20) is 
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(3.4.22) 



-V 4<^vT| 



The second term of F 2 (s) is evaluated by convolution. Except 
for a constan 



for a constant, it is the first term of f ^ ( t ) convolved with 



or 






(3.4.23) 



The expression (3.4.23) can be evaluated as 



\wUi b \_ 



_ (oi. - Js — Cc^(?J^-v)tl 

— ^ (c^y- v - 1 



o 0 -\ 



(3.4.24) 



where 



w 

k 



\s - 









Q, 


4 


4 




W 


^'V 5 

(o 


(^V 0* 

v 4 


4 




(O b W> 




C^o-0 3 


k> 


Vb 


4 


_ , 4 






C^VD"- 





^6 



r \ 



i^X \ 






(3.4.24a) 
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Except for a constant, the last term of F 2 (s) is 
simply the term of (3.4.24) convolved again with 1_ -V 

The result is found in Appendix C. 

Clearly the method is too cumbersome to carry beyond 
this second approximation. The approximate solution by the 
Brady-Baycura technique is 






(3.4.26) 



3 . Comparing the Solutions 

The time plots of Figures 16, 17, 18 and 19 exhibit 
the relative accuracy of the approximations to the Runge- 
Kutta solution, for various combinations of x q , B and 

From Table IV it is apparent that the approximations 
are good only for the first initial time units. Pipes method 
gives only a slightly better approximation than the Brady- 
Baycura technique. For the case (S q = 5.0, it is observed 
that the approximations by both methods improve. This is 
in keeping with the original restriction of the system to 
forcing functions far from resonance. 

E. NONLINEAR SPRING WITH CUBIC TERM AND FORCING FUNCTION 
Assume the describing equation to be 

tT* + X VX 3 -- 

^ (3.5.1) 



with initial conditions 
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*<£) 




FIGURE 16. 



Time solution of x + x + x 

x = 0.5 
o 



0.5 cos 1 . 5t , 
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* 




FIGURE 17. 



Time solution of x + x + x 
x = 0.5 



0 . 5 cos 5 . Ot , 
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x(i) 




FIGURE 18. 



Solution to x + x 




0.5 cos 1.5t, x Q 0.25 
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x<§ 




FIGURE 19. 



Time solution of x + x + x 

x = 0.5 
o 



0.25 cos 1 . 5t , 



60 



TABLE IV 





X 

0 

II 

w 

II 

o 

• 

<J1 


o Z1 ' 5 






Time 


Runqe-Kutta 


Pipes 




Brady 


0.000 


0.50000 


0.50000 




0.50000 


0.050 


0.49969 


0.49969 




0.50015 


0.100 


0.49875 


0.49875 




0.50053 


0.250 


0.49208 


0.49205 




0.50217 


0.500 


0.46719 


0.46657 




0.49946 


1.000 


0. 35328 


0.34448 




0.41322 


1.500 


0. 13099 


0.09665 




0.14546 




x “ B ~ 0.5, 
o 


- 5.0 
o 






0.000 


0.50000 


0.50000 




0.50000 


0.050 


0.49968 


0.49968 




0.49991 


0. 100 


0.49870 


0.49870 




0.49955 


0. 250 


0.49034 


0.49030 




0.49466 


0.500 


0.44401 


0.44340 




0.45529 


1.000 


0. 17816 


0.17054 




0.17995 


1.500 


-0.12313 


-0.14487 




-0.18418 




x = 0. 25 , B 
o 


II 

o 

• 

U"i 

0 

II 


1.5 




0.000 


0.25000 


0. 25000 




0. 25000 


0.050 


0.25023 


0. 25023 




0. 25041 


0.100 


0. 25093 


0.25093 




0.25161 


0.250 


0. 25563 


0. 25563 




0. 25960 


0.500 


0. 26987 


0.26979 




0.28343 


1.000 


0.29115 


0. 28972 




0.32447 


1.500 


0.22963 


0. 22239 




0. 26635 




x q = 0.5, B = 


0.25, Q = 


1.5 




0.000 


0.50000 


0.50000 




0.50000 


0.050 


0.49937 


0.49938 




0.49972 


0.100 


0.49750 


0.49750 




0.49881 


0.150 


0.49438 


0.49438 




0.49722 


0.250 


0.48444 


0.48441 




0.49164 


0.500 


0.43861 


0.43801 




0.46008 


1.000 


0.26682 


0. 25904 




0. 29380 


1.500 


0.015343 


-0. 01124 




-0.01725 
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XCoV x e 

kC°\ = o 



(3.5.1a) 



1 . 



a . s . 

i 



Solution by Pipes Method 

The solution proceeds in the usual way be identifying 

C\ = T 5 >\ 

t 

^ ~ ' ( (3.5.2) 

W = l 
0(0 - 



Using equation (2.1.7), the first approximation A^, is found 

/ 

by Laplace transform techniques. 









(3.5.3) 



A v co - 












(3.5.4) 



The solution is the same as for the quadratic case 



A,C0- fe-,c- o ^sL - 



(3.5.5) 



where 



k ,= Xc -V 

V - ^ 






(3.5.6) 
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The coefficient A 2 (t) was seen previously to be zero. 
Solving for (t) is straightforward and the result is 

A 3 uv- + 

L v - n 



- -3^ - 3 w,\C: . wt \ , » 

4\CZ^-'T-D 



-V 3^\ co^to.v^ 

3^1^ {Lo^Cg-^t^-f _ -3 C-osfz-^V'H 

' 4 t Zu v b-Q 0 

4- V> 4 “] 

\ 1 Cos>2.coji 

4^^-0 J 



where 

W = SV? v > k.wt ■) 

^ *" ( 

(3.5.8) 

W A - ^ 4 3/, 

The solution of the undriven system of Chapter III, Section 
C should demonstrate the unfeasibility of Carrying the 
approximation beyond A^(t). Furthermore, the first secular 
term in A^(t) will cause the solution to diverge. 

2. Solution by Brady-Baycura Technique 

Taking Laplace transforms of equation (3.5.1) and 
impressing the initial conditions, the system equation 
becomes 

^ XV) -V 0?* b)Xts) - *SX 0 -v — (3.5.9) 

S*-V uV- 

o 
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Using the algebraic reversion of series technique to solve 
for X(s), let the left side of equation (3.5.9) be a series 



-y \^X O) A k 3 X Cs) \ 



(3.5.10) 



then 



XfcV 






-v . . . 



(3.5.11) 






o. 






\ it-v-cCy*- 



V * . , (3.5.12) 



Where from Ref. 9 






1_ _ 



i* -V\ 



(3.5.13) 



^ - - 






= D 



(3.5.14) 






- 



k>A (yvO A 



(3.5.15) 



Then to first order approximation 



xuv- 






(3.5.16) 









Jb -v \ 






(3.5.17) 



and 






(3.5.18) 



64 



where 



Jk_ 

u^-y 



(3.5.18a) 



k = ^ 

o 



w 

* z — 

U^-_\ 



As expected, the result is equivalent to Pipes method for 
the first approximation. 

Since B 2 is zero, the third term, F 3 (s), of the 
series is found 









V 









(3.5.19) 












S — 4. - o 

t^'T <£*■'?(£<"%) I) 



-t 






texuaSKs 



V X A 



(3.5.20) 



The inverse transform of the first term as given in 
equation (3.3.18). 







(3.5.21) 



Except for the constants, the successive terms of F 3 (s) must 

be evaluated by successive convolution of — sin^\ t with 

^o 0 

the right side of equation (3.5.21). While continuing the 
evaluation of F 3 (s) will lead to terms of cos t, modulation 
terms and higher resonance terms, it will also contribute 
numerous secular terms, and thus cause the solution to di- 
verge more rapidly. Again the technique has lost its desir- 
ability by providing cumbersome transforms. 
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Using the first term only of F^s), the solution is 

given by 

xCAVs.k.Cod - x? fis-Ui-A 

4 «& \_ 



+ ?> a\, (3.5.22) 

— t. -b\rrt 1+ . . , ' 



3. Comparing the Solutions 

Figures 20, 21, 22 and 23 are graphs of the approxi- 
mate solutions and the true solution for various combinations 

of x , B anduA . Table V lists selected values from the 
o o 

solutions . 

The graphs demonstrate that except in the extreme 
case of Figure 20 where x q = B and <£) = 1.5, that Pipes 

Method provides only a slightly better approximation than 
the Brady-Baycura Method. This is due to the number of 
terms generated by the Pipes Method. Figure 21 shows that 
the condition of the system operating far from resonance. 
Pipes Method gives a better solution for a longer period of 
time. Both approximations improve as x q and B become small. 
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x<&) 




FIGURE 20. 



Time solution of x + x + x 
x = 0.5 



0.5 cos 1.5t, 
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FIGURE 21. Time solution of x + x + x 0.5 cos 

x^ = 0.5 , 

o 



68 






*■ 





< 3-00 1-00 



ISO z.°c> 



FIGURE 22. 



Time solution to x + x + x 
x =0.25 



0.5 cos 1.5t, 
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XCO> 



— \ 




FIGURE 23. 



Time solution to x + x + x 

x = 0 . 5 
o 



0.25 cos 1.5t, 
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TABLE V 





x - B ~ 0.5, 
o 


= . 1.5 






Time 


Runqe-Kutta 


Pipes 




Brady 


0.000 


0.50000 


0.50000 




0.50000 


0.050 


0.49984 


0.50017 




0.49984 


0. 100 


0.49937 


0.50068 




0.49937 


0.150 


0.49857 


0.50150 




0.49858 


0.250 


0.49594 


0.50386 




0.49599 


0.500 


0.48209 


0.51016 




0.48280 


1.000 


0.40483 


0.47304 




0.41464 


1.500 


0.22075 


0. 28681 




0. 26156 




X Q ~ B = 0.50, 


= 5.0 






0.000 


0.50000 


0.50000 




0.50000 


0. 050 


0.49984 


0.49988 




0.49984 


0.100 


0.49932 


0.49949 




0.49933 


0.150 


0.49834 


0.49870 




0.49835 


0. 250 


0.49420 


0.49516 




0.49425 


0.500 


0.45886 


0.46195 




0.45941 


1.000 


0.22743 


0. 23117 




0.23216 


1.500 


-0.04477 


-0.04985 




-0.02377 




x q = 0. 25 , B = 


0.50, 


= 1.5 




0.000 


0. 25000 


0.25000 




0. 25000 


0.050 


0. 25029 


0. 25034 




0. 25029 


0.100 


0. 25116 


0. 25135 




0. 25117 


0.250 


0. 25709 


0.25816 




0. 25711 


0.500 


0. 27572 


0. 27884 




0. 27596 


1.000 


0. 31407 


0.31302 




0. 31748 


1.500 


0. 27665 


0.25064 




0. 29059 




x = 0.50, B = 
o 


0. 25 , 


= 1.5 




0. 000 


0.50000 


0.50000 




0.50000 


0.050 


0.49953 


0.49962 




0.49953 


0.100 


0.49812 


0.49846 




0.49813 


0. 150 


0.49578 


0.49653 




0.49579 


0.250 


0.48829 


0.49029 




0.48831 


0.500 


0.45336 


0.45969 




0.45362 


1.000 


0.31626 


0. 32228 




0.32072 


1.500 


0.09741 


0.07727 




0. 12178 
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IV. CONCLUSIONS . 



It has been demonstrated that these various methods of 
obtaining approximate solutions to nonlinear systems with 
initial conditions are accurate only during the early time 
intervals. Further, it was shown that this degree of accuracy 
may be satisfactory for engineering applications. 

It has also been shown that the Brady-Baycura technique 
of nonlinear transforms usually provides a more accurate 
solution; however, the ease and facility of Pipes Method 
might generally make this a more desirable method. 

The method of handling so-called secular terms is still 
unresolved. .It is suggested that they may be eliminated by 
letting the frequency also be a series such that the con- 
stants of the series are chosen to eliminate resonance con- 
ditions and thus eliminate the secular terms. The method 
would be similar to that given by Pipes in Ref. 12. 

Both methods show that the higher order harmonic terms 
exist in' the solution; however, the subharmonic terms fail 
to occur in the first few terms of the approximations. 

Since the problem is an initial value problem, it appears 
that the solution could be improved by stopping the problem 
solution before it diverges from the true solution (t .5), 
then restarting with a new set of initial conditions. The 
effects of x(o) would have to be investigated. The value of 
such a method seems dubious since Runge-Kutta analysis on the 
digital computer seems an easier method. 
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APPENDIX A 



INVERSE LAPLACE TRANSFORM OF 

(s 2 + 1) 



The derivative technique is used : 
From Ref. 10, 
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o = 






(A- 4) 



4- V* 

where n is the n derivative. Thencombining equations 
(A-4) and (A-3) 




(A-5) 



(A- 6) 
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APPENDIX B 



EVALUATION OF A 3 IN FORCED QUADRATIC SYSTEM 



From equation (2.1.5) A^ becomes 

Ajiv - 



(B-l ) 



After substituting for a^, A^ and A 2 and combining terms, 

- W^sc^-t 

t> u -v\ 6 (o' 

+- -v v WJk + Cx>^(j^Y>’H 



(B-2) 



where 
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Except for the first term of A^ , the transforms are all of 
the form : 

i.-nA 

(B-4) 



\) (sW) 



For which the inverse transform is known: 
J \ tubs 



LO'-\ 



(B-5 ) 



The first term has the transform: 

which has the inverse 

£ 

Using expressions (B-5) and (B-7) , A^ becomes 
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APPENDIX C 



SOLUTION OF LAST TERM OF F^s) BY BRADY-BAYCURA METHOD 

Let the expression in parenthesis of (3.4.24) equal Z^(t) 

and sin t equal Z_(t), then 
U> 0 ° 2 V 

^ ■£.(& V-f Jl-0 + . - . (C-l) 

thus 

i£_ . .. 

3*.0)‘L J (C-2) 

or six convolutions make up the last term. Since each of 
these convolutions is of the type previously evaluated, the 
solution is now straightforward and the last term becomes 
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